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(54) Improved method and apparatus for determining dip angle, and horizontal and vertical 
conductivities 

(57) The present invention is directed to a method 
and apparatus for determining the dip angle and hori- 
zontal and vertical conductivities of a wellbore during 
drilfing operations- If the dip angle is not known, the 
present invention allows for the use of muftffrequency 
measurements and a model which maps dielectric per- 
mitivities to conductivities to calculate the horizontal 
conductivity, the vertical conductivity, the horizontal die- 
lectric permitivity, the dielectric perrnrtivrty. and the dip 
angle. If the dip angle is known, the present invention 
allows a muttifrequency tool to be utilized alone to cal- 
culate the horizontal conductivity, the vertical conductiv- 
ity, the horizontal dielectric perrnitivrty, and the vertical 
dielectric permitivity. If the dip angle is not known, a sin- 
gle frequency resistivity tool (as opposed to a rnuttifre- 
quency tool) can be utilized in corftwiation with a model 
which maps dielectric permitivities to conductivities to 
calculate the horizontal conductivity and the horizontal 
dielectric permitivity. 



S2.S4 



USEEQUNlONS/S'ftCFCfWMO 



usE^to.aaaWGeeanDjofM 



HO 



ODSEBCUGH 



to 



HO. 2 



WmKJby X««* (UK) Business Services 
2.16.3/34 



EP 0840 142 A2 



Descripti n 

BACKGROUND OF THE INVENTION 
5 1. Field of the Invention: 

The present invention relates in general to well-logging operations, and is particularly useful in anisotropic forma- 
tions. 

io 2. Description of the Prior Art: 

In oil and gas exploration, measurement-while-drilling or wireline systems are utilized to develop logs of the earth 
formations surrounding the borehole. In general, logging packages include: a transmission system that energizes the 
earth formations with either an electromagnetic field or currents, and at least one receiver system for monitoring the 
is effect that the earth formations (and borehole) have on the field or current 

The electrical properties of the earth formations provide information about the geologic materials that make up the 
formations, and about their likely oil. gas. and water content The most useful electrical attribute is resistivity (or con- 
ductivity). However, the dielectric property of the earth formations is also of interest 

Some earth formations are isotropic, which means that the earth formation does not preferentially conduct electric- 
20 ity in any particular direction. Other formations are anisotropic, which means that the earth formation may have a 
greater conductivity in one particular direction than in another particular direction. This is particularly true in many sed- 
imentary geologic formations in which current flows more easily in the direction parallel to the bedding planes than 
transversely to them due to the fact that a number of mineral crystals are flat or elongated in shape, and they naturally 
took an orientation parallel to the plane of sedimentation when they were deposited. 
25 In oil and gas exploration, the anisotropy of earth formations can be expressed as a horizontal conductivity, and a 
vertical conductivity which is transverse. Alternatively, an anisotropy coefficient (X) can be utilized to quantify the anisot- 
ropy in accordance with Equation No 1 : 



The determination of horizontal and vertical conductivities is complicated by directional drilling which generates 
boreholes that are at an angle relative to the earth formations. This angle is called relative dip angle and is often not 
known precisely. 

40 Others have proposed techniques for determining conductivity in anisotropic formations. In an article by T. Hagi- 
wara entitled "A New Method to Determine Horizontal-Resistivity in Anisotropic Formations Without Prior 
Knowledge of Relative Dip", presented on June 1 6-1 9. 1996. at the SPWLA 37th Annual Logging Symposium, a tech- 
nique is proposed which allows for the calculation of horizontal conductivity without any prior knowledge of the dip 
angle. This may reduce the dependence of logging-while-drilling operations on survey instruments. 

45 One problem with the approach suggested by Hagiwara is that influence of the dielectric properties of the formation 
is largely ignored. In relatively high resistivity geologic formations (like oil-bearing sands) the impact of the dielectric 
effect increases In influence on the measurements from which conductivity is calculated. Also, the dielectric effect has 
a greater impact for high frequency interrogation signals. 

50 SUMMARY OF THE INVENTION 

It is one objective of the present invention to allow for the simultaneous measurement of horizontal and vertical 
resistivities, horizontal and vertical dielectric permitivity (or "constant"), and the dip anc^e using multiple frequency 
resistivity tools and using anisotropic formation models that relate dielectric permitivity to conductivities. 
55 Alternatively, if the dip angle is known, the conductivities and dielectric constants can be measured using muttifre- 
quency tools without any anisotropic formation model. The measurements can also be used to verify the validity of the 
anisotropic formation models which relate horizontal and vertical formation conductivities to the dielectric permitivity. 
This is not practical with single frequency tools. 
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EQUATION NO- 1: 




35 
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Alternatively, if the model which relates horizontal and vertical conductivities t the dielectric effect is known, this 
invention also allows one to obtain the horizontal resistivity, as well as the horizontal dielectric permitivrty, using single 
frequency tools without knowing the dip angl in high resistivity zones or with high frequency interrogation signals. If the 
dip angle is known, the vertical resistivity and dielectric permitivity can also be measured. 

5 

BRIEF DESCRIPTION OF THE DRAWINGS 

The novel features believed characteristic of the invention are set forth in the appended claims. The invention itself, 
however, as well as a preferred mode of use, further objectives and advantages thereof, will best be understood by ref- 
it? erence to the following detailed description of an illustrative embodiment when read in conjunction with the accompa- 
nying drawings, wherein: 

Figure 1 depicts a drillstring which includes a logging-while-driHing system. 
Figure 2 is a flow chart representation of the implementation of the present invention. 
is Figures 3 through 7C depict a logging^iledrilltng tool that can be utilized to implement the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

It has been shown by Hagiwara that using resistivity tools with axial<fipole antennas, one can measure horizontal 

20 resistivity and a function of vertical resistivity and relative dip angle. Without additional information, the vertical resistiv- 
ity and relative angle cannot be uniquely determined, no matter how many antenna soarings are used. 

However, this is only true for low frequency tods. At higher frequencies, where the dielectric effect becomes impor- 
tant, even the horizontal resistivity cannot be determined with a single frequency measurement in deviated wells. The 
antenna readings are affected by five formation parameters: two resistivities, two dielectric constants (or "perrnrtivities"), 

25 and dip angle. Four independent parameters derived from the five formation parameters completely determine the 
antenna readings. They are horizontal resistivity, horizontal dielectric permitivity. a function of vertical resistivity and rel- 
ative dip angle, and the same function of vertical dielectric permitivity and relative dip angle. With a single frequency, 
the contributions from displacement and conduction currents in the formation to the measurements cannot be, as a 
practical matter, separated. Since the effects of dielectric permitivity and resistivity on wave number vary differently with 

30 frequency, in accordance with the present invention, with multiple frequency measurements one can obtain the four 
parameters. However, these four atone cannot determine all five formation parameters. 

One other function or equation is necessary in order to determine all five parameters. The function or equation that 
allows the determination of all five formation parameters is the relationship between conductivity and the dielectric 
effect. The factors that cause macroscopic conductivity anisotropy necessarily cause dielectric anisotropy. The aniso- 

35 tropic dielectric permitrvities and conductivities are not independent of each other. Once the relationship between them 
is established for a given formation type, one can use the four measurements, plus a formation anisotropy model (or 
chart, or formula), to uniquely determine all five parameters. Of course, if the relative dip angle is known, the resistivities 
and dielectric permitrvities can be determined without any formation models. 

As depicted in Figure 1 , in accordance with the preferred embodiment of the present invention, a logging-while-drill- 

40 ing system 11 is provided and coupled to drillstring 21 which extends into borehole 23. As is shown, borehole 23 
extends at an angle through earth formation 25 which, for example, is composed of alternating layers of sand and shale 
deposits. It includes at least one transmitter 1 3. 1 5 for generating interrogating electromagnetic fields. Preferably, the at 
least one transmitter 1 3, 1 5 can generate interrogating electromagnetic fields having different frequencies. The logging- 
while-drilling system 1 1 also includes at least one receiver for detecting the effect of the earth formations (and borehole) 

45 on the interrogating electromagnetic fields. Preferably, a differential receiver is utilized which includes two spaced-apart 
receiver antennas 17, 19 and reception circuitry for detecting current induced in the receiver antennas 17, 19. 

As is conventional, logging-while-drilling system 11 is capable of measuring any of the following electrical attributes 
which can be utilized to measure any combination of the following, utilizing at least two interrogation frequencies: 

so (1) the amplitude attenuation of the interrogating electromagnetic field; 

(2) the phase shift of the interrogating electromagnetic field; and 

(3) the in-phase and out-of -phase components under phase quadrature analysis. 

In addition to the measurements, the present invention requires that a model (formula, chart table, eta) be devel- 
65 oped that relates (or "maps") the dielectric "constants" or "permitivities" of a particular type of earth formation to partic- 
- -ular anisotropic conductivities. ~ r 

The present invention allows for the following particular uses: 



3 



EP 0840 142 A2 



1. If the dip angle is not known, the present invention allows for the use of (1) the multifrequency measurements, 
and (2) the model which maps dielectric perrnrtivities to conductivities, to calculate: 

(1) the horizontal conductivity; 
5 (2) the vertical conductivity; 

(3) the horizontal dielectric permitivrty; 

(4) the vertical dielectric perrnrtivity; and 

(5) the dip angle. 

io 2. H the dip angle is known, a multifrequency resistivity tool can be utilized alone (that is, without the availability of 
the model which maps dielectric perrnrtivities to conductivities) to calculate: 

(1) the horizontal conductivity; 

(2) the vertical conductivity; 

is (3) the horizontal dielectric permrtivity; and 

(4) the vertical dielectric permrtivity. 

3. H the dip angle is not known, a single frequency resistivity tool (as opposed to a multifrequency tool) can be uti- 
lized in combination with the model which maps dielectric per mativi ties to conductivities to calculate: 

20 

(1) the horizontal conductivity; and 

(2) the horizontal dielectric permativity. 

How these commercially advantageous uses can be obtained utilizing the present invention will now be described. 
25 The tog response of an induction tool with axial dipole antennas to an anisotropic formation is described by Equa- 
tion No. 2: 

EQUATION NO. 2: 

30 i 



[-2(1-1. k k . ^ +<K ■ L) (e« s - e"*0] 



35 wherein the induced voltage V is expressed in terms of three parameters of Equations 3, 4. and 5. as follows: 

P EQUATION NO. 3: 



k,S + k h Lfi 



EQUATION NO. 4: 



cos 2 6 • L 



so EQUATION NO. 5: 



I P 

= Jcos^+^Y-sin^ 



wherein: 
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0 represents the dip angle 

k|! represents the complex wave number in th horizontal direction 

represents the complex wave number in the vertical direction 
L represents the spacing between a transmitter and a receiver 
5 S represents the rescaled distance in an anisotropic medium 

The values for kj, and can be determined from the following Equations No. 6 and 7: 
EQUATION NO. 6: 



10 



15 EQUATION NO. 7: 



20 



25 



30 



Equation No. 5 can be restated as Equation No. 8: 
EQUATION NO. 8: 



P 2 ( 



If the relative magnetic permeability is assumed to be 1 (which is a fair assumption for non-ferrous formations), 
Equation No. 8 can be restated as Equation No. 9: 

35 

EQUATION NO. 9: 



In formations where the CRIM model (which is a model that relates resistivity to dielectric permitivity, as described 
45 in US. Patent 5,144,245, 'Method for Evaluating a Borehole Formation Based on a Formation Resistivity Log Gener- 
ated by a Wave Propagation Formation Evaluation Tool", to M. Wisler, and owned by Baker Hughes Incorporated, and 
incorporated herein by reference) applies, the dielectric permitivity is related to the formation conductivity through 
Equation No. 10, as follows: 

50 EQUATION NO. 10: 



55 where e is the dielectric constant a is the formation conductivity, c and a are constants. At 2 MHz, c ■ 21 0 and 

a «= -.43. At 400 KHz, a is slightly Afferent and c is about twice of its value at 2MHz: Since the dielectric effect at 400> ^ - 
MHz is much smaller than that of 2 MHz, it is a good approximation to assume both c and a as being independent of 
frequency for MPR tool measurements. This approximation is made for simplicity of arguments. Any known frequency 



5 



EP 0 840 142 A2 

dependency used in the CRIM model would only increase the mathematical complexity of the application of the inven- 
tion. It would not change the essence of the invention. In the present invention th "model" which releases relates resis- 
tivity to dielectric permitivity may be a data file maintained in electronic or semiconductor memory which is accessible 
through machine readable program instruction, or it may constitute one or more equations and/or constants also stored 
s in memory. 

The compensated long spacing phase differences and amplitude ratios at 400 KHz and 2 MHz are p400k a400k, 
p2m, a2m. They can be used through an inversion method to obtain the four parameters in accordance with Equations 
Nos. 11-14. as follows: 

10 EQUATION NO. 11: 

ml=e A 



EQUATION NO. 12: 



EQUATION NO. 13: 



m2 = a h 



nt3 = (e v -e M )sm 2 Q 



EQUATION NO. 14: 

30 /«4 = (a„ -cr A )sin 2 9 

35 The inversion involves four parameters (and is thus "four-dimensional") and can be extremely complicated. An 
effective inversion method is an iterative one which reduces the four-dimensional inversion into a series of two-dimen- 
sional inversions which have been used in the prior art Since the dielectric effect on the tower frequency measurements 
is small, one starts by assuming that there is no dielectric effect on the 400 KHz measurements. The two lower fre- 
quency measurements are used to obtain m2 and m4. Then the two 2 MHz measurements are used to determine the 

40 remaining two parameters, ml and m3. The non-zero ml and m3 are then used to refine m2 and m4 by a new 400 KHz 
inversion. The newly obtained m2 and m4 are used in the 2 MHz inversion to update ml and m3. This process contin- 
ues until a predetermined accuracy is achieved. The attached flow chart of Figure 2 shows how the 4-d inversion proc- 
ess works. 

After the inversion, we have the following Equations Nos. 15-18. 

45 



50 



65 
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EQUATION NO. 15: 



EQUATION NO. 16: 



10 



ml=e A 



m2 = a A 



EQUATION NO. 17: 

is m3 = (e„ - ejsin 2 8 



EQUATION NO. 18: 



20 



m4 = (a„ -a A )sin 2 9 



where ml. m2, m3, m4 are the results of inversion. The horizontal dielectric perrnitivity and resistivity are determined. 
Using the CRIM model we have Equation Nos. 19 and 20. 

25 

EQUATION NO. 19: 

m3 = c(cr; -a;)sin 2 G = c(a tf a - m2 a )sva 2 Q 

30 

EQUATION NO. 20: 

m4 = (<j v -ajsin 2 9 =(o v - ml) sin 2 Q 

35 

The ratio of the above two equations gives us Equation No. 21. as follows: 

40 



EQUATION NO. 21: 



m3 <jc*-m2 a ) 
m4 o-m2 



45 



This equation is used to calculate o v through a trivial one-dimensional inversion. Using the CRIM model. is cal- 
culated. So we have determined e^e^o^o^. Using either Equation 17 or 18. the dip angle is determined. 
so The apparatus that can be utilized to implement the present invention will now be described with reference to Fig- 
ures 3 through 7C. 

TRANSMISSION AND RECEPTION SYST1EMS; Figure 3 provides a block diagram view of an exemplary logging 
tod 201 constructed in accordance with the present invention. Logging tool 201 includes upper transmfters 203. 205. 
lower transmitters 207, 209. and intermediate series resonant receiving antennas 211. 213. Central processor 215 is 
«5 preferably a mkxoprocessor device which is utilized to coordinate the operation of the components of logging tool 201 . 
t record and process the data obtained from measurements made by intermediate series resonant receiving antennas 
211. 213, and to interact with the mud pulse telemetry data transmission system carried in the adjoining drill collar 
member. Processor 217 is provided and dedicated for the control of numerically controlled oscillator 223. Processor 
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21 9 is provided and dedicated for the control of numerically controlled oscillator 225. Central processor 21 5 communi- 
cates with processors 217, 219 via data buses 241 , 243 respectively. Numerically controlled oscillators 223, 225 are 
adapted to receive a binary command signal as an input and to produc an analog output having particular frequency, 
phase, and amplitude attributes. Th frequency, phase, and ampTrtude attributes are determined at least in part by the 

5 command signals applied from processor 217. 219 1 the input of numericaily controlled oscillators 223. 225. and the 
data contained in various registers within numerically controlled oscillators 223. 225. Numericaily controlled oscillators 
223, 225 provide the analog signal to transmitting circuits 227, 229 respectively. 

Receiving antennas 211. 213, communicate through analog receiving circuit 231 with the first and second data 
input channels of a digital signal processor 221. The digital signal processor 221 receives data at the first and second 

10 inputs after it is converted from analog form to digital form by analog-to-digital converters 220, 222, and records the 
data elements In a circular memory buffer. Central processor 21 5 pulls data from the buffers in a prescribed and prede- 
termined manner in order to sample the current which is generated in receiving antennas 211, 213 in response to the 
propagation of electromagnetic signal through the adjoining formation. As is conventional, the resistivity of the forma- 
tion surrounding the logging tool 201 may be determined by either (1) determining the amplitude attenuation of an elec- 

15 trorrtagnetic wave propagating through the formation adjoining receiving antenna 21 1 and receiving antenna 213, or (2) 
by determining the phase shift between the electromagnetic signal propagating through the formation adjoining receiv- 
ing antenna 21 1 and 213, or from both. These measurements comprise a relative measurement of the amplitude atten- 
uation and a relative measure of the phase shift 

The present invention also allows other techniques for quantifying the electromagnetic field which propagates 

20 through the formation surrounding logging tool 201 . Since precise control can be obtained with the present invention 
over the frequency, phase, and amplitude of the electromagnetic wave generated by transmitting antennas 203, 205, 
207, and 209, the present invention allows the measurement of the absolute amplitude attenuation of electromagnetic 
signal between any particular transmitting antenna 203, 205, 207, and 209 and any particular receiving antenna 21 1 , 
213. Furthermore, the logging tod 201 of the present invention allows for the absolute measurement of the phase shift 

25 of an electromagnetic signal between any particular transmitting antenna 203, 205, 207. 209 and any particular receiv- 
ing antenna 211, 213. Prior art devices do not allow such optional techniques for determining amplitude attenuation and 
phase shift since prior art devices are unable to determine easily and precisely the frequency, phase, and amplitude of 
a signal generated at any particular transmitting antenna 

The operation of numerically controlled oscillators 223. 225 is clocked by the output of reference clock 237, which 

ao is preferably 1 2 megaHertz. The operation of receiving circuit 231 is controlled by the output of numerically controlled 
oscillator 231 . which is also clocked by the output of the reference clock 237, which is 12 megaHertz. Thus, a clocking 
pulse is provided to numerically controlled oscillator 223. 225 at a frequency identical to that which is provided to 
numerically controlled oscillator 223, which establishes the operating frequency of receiving circuit 231. Digital signal 
processor 221 is clocked by the output of divide-by circuit 239. and thus samples the output of receiving circuit 231 at 

35 a particular frequency which is much less than that utilized to energize transmitter antennas 203, 205, 207 and 209. 

Numerically controlled oscillator 233 produces a phase-locked sine-wave signal with a center frequency of 1 .995 
Mhz, that is used as a local oscillator signal by receiving circuit located in electronics housing 73. 

Reference is now made to Figure 4. The overall function of the circuitry depicted in block diagram and schematic 
form in Figure 4 is to respond to the focal oscillator signal and one of the two receiver coil output signals to produce a 

40 receiver phase output signal relative to the transmitter and a receiver amplitude output signal. A conventional pre-amp 
circuit generally indicated at 271 responds to the receiver pick-up signal and its output is applied to a mixer circuit 
arrangement generally indicated at 273. Mixer circuit arrangement 273 includes an integrated circuit 275 that suitably 
is implemented by an integrated circuit manufactured and sold by Motorola and other companies under the designation 
MC1596. 

45 Because the frequency of the pick-up signal and the local oscillator signals are phase-locked to a common fre- 
quency reference and differ by 6 Khz, the intermediate frequency (IF) produced by mixer circuit arrangement 273 is at 
6 Khz. A band pass tuning circuit arrangement generally indicated at 277 passes the 6 Khz IF signal to an amplifier cir- 
cuit arrangement generally indicated at 279. An active band pass fitter circuit arrangement generally incScated at 281 
provides further band pass filtering and provides a signal to an analog-to-digital converter, which supplies a digital input 

50 to a particular input channel of digital signal processor 221 (of Figure 3). 

Figures 5A and 5B is a block diagram view of the numerically-controlled oscillators 223. 225, 233 of Figure 3. 
Since the numerically-controlled oscillators are identical, only numerically-controlled oscillator 223 will be discussed 
and described. In the preferred embodiment of the present invention, numerically-controlled oscillator 223 comprises a 
CMOS, DOS modulator manufactured by Analog Devices of Norwood. Massachusetts, which is identified by Model No. 

65 AD7008. The numerically-controlled oscillator 223 includes a thirty-two bit phase accumulator 301 . a sine and cosine 

- look-up table 303, and a ten-bit digital to analog converter 305. Clock input307*isprovided to receive a docking 6ignal 

from a device which is external to the numerically-controlled oscillator 223. The particular numerically-controlled oscil- 
lator of the present invention is adapted to accept clock rates as high as twenty megaHertz to fifty megaHertz, but can 
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accommodate much lower clock rates. Th device purports t have a frequency accuracy which can be controlled to 
one part in four billion. Numerically-controlled oscillator 223 includes a thirty-two bit serial register 309 which receives 
serial data at serial data input pin 31 1 . which is clocked into th register in accordance with a dock signal which is sup- 
plied to serial dock input 313. A thirty-two bit parallel register 313 is also provided which receives parallel binary data 

5 from MPU interface 315. Data bus 317 includes sixteen digital input pins identified as DO through D15. The chip select 
pin 321 is utilized when writing to the parallel register 313. The write pin 319 is also utilized when writing to the parallel 
register 309. The transfer control address bus 323 is utilized to determine the source and destination registers that are 
used during a transfer. A source register can be either the parallel assembly register 313 or the serial assembly register 
309. The destination register can be any one of the following registers: the command register 325. the FREQO register 

10 327, the FREQ1 register 329, the phase register 331 , the IQMOD register 333. The command register is written to only 
through the parallel assembly register 313. The contents of the command register determine the operating state of the 
nurnerically-controOed oscillator 223. In the preferred device utilized in the present invention, the command register is 
a four bit register. The content of this register determines the operating state of the numerically-controlled oscillator. 
During logging operations, the logging apparatus of the present invention is programmed to provide commands from 

15 processors 215, 217, 219 (of Figure 3) with eight-bit commands, so the "CW bit is 0. Normal operation is desired, so 
the "CR1 " bit is 0. In the present invention, amplitude modulation is bypassed, so the "CR2" bit is 0. In the present inven- 
tion, the synchronizer logic is enabled so the "CR3" bit is 0. The FREQO register 327 defines the output frequency of 
the numerically-controlled oscillator 223. when the FSELECT pin is 1 , as a fraction of the frequency of the clock signal 
applied to clock pin 307. The FREQ1 register 329 defines the output frequency of the numerically-controlled oscillator 

20 223, when FSELECT equals 1 , as a frequency of the dock signal applied to clock pin 307. The contents of the phase 
register 331 are added to the output of the phase accumulator 301. The IQMOD register 333 is not utilized in the 
present invention. 

The operations can be performed with the registers by supplying command signals to transfer control address bus 
323. Three basic operations can be performed. The contents of the parallel assembly register 313 can be transferred 

25 to command register 325; the contents of the parallel assembly register can be transferred to a selected destination reg- 
ister; and the contents of the serial assembly register 309 can be transferred to a selected destination register. 

The load register pin 335 is utilized in conjunction with the transfer control address bus 323 to control loading of 
internal registers from either the parallel or serial assembly registers 309, 313. The test pin 337 is utilized only for fac- 
tory testing. The reset pin 339 is utilized to reset the registers. The reset pin in particular is utilized to clear the command 

30 register 325 and all the modulation registers to 0. The current output pins 341 , 343 are utilized to supply an alternating 
current to a selected end device. In the particular embodiment of the present invention, only one of these outputs is uti- 
lized for a particular transmitting antenna, since one current is the cornpliment of the other current The compensation 
pin 342 is utilized to compensate for the internal reference amplifier. The voltage reference pin 343 can be utilized to 
verride an internal voltage reference, if required. The full-scale adjust pin 345 determines the magnitude of the full 

35 scale current at output pins 341 , 343. The ground pin 347 provides a ground reference, while the positive power supply 
pin provides power for the analog components within numerically-controlled oscillator 323. The frequency select pin 
351 controls frequency registers FREQO register 327 and FREQ1 register 329, by determining which register is used 
in the phase accumulator 301 by controlling multiplexer 353. Trie contents of phase register 331 is added to the output 
of phase accumulator 301 at sumer 355. The IQMOD registers 333 are provided to allow for either quadrature ampli- 

40 tude modulation or amplitude modulation, so the sine and cosine outputs of look-up table 303 are added together at 
sumer 357, and are unaffected by the IQMOD registers 333. The output of sumer 357 is provided to digital-to-analog 
converter 305, which creates an analog signal having a frequency which corresponds to either the contents of the 
FREQO register 327 or the FREQ1 register 329. a phase which is determined by the output of sumer 355 which is pro- 
vided as an input to look-up table 303. and an amplitude which is determined by full scale control 359 which is set by 

45 full scale aojust pin 345 and reference voltage pin 343. Therefore, the numerically-controlled oscillator of Figures 5A 
and 5B can provide an analog output having a precise frequency attribute, phase attrfoute, and amplitude attribute. 
Since the device is extremely accurate, it is possible to provide a driving current for the transmitting antennas 203, 205, 
207, 209 of Figure 3 which is controlled precisely. In the preferred embodiment of the present invention, one of trans- 
mitting antennas 203, 205 is operated at 400 kiloHertz, while the other of transmitting antennas 203. 205 is operated at 

so 2 mega Hertz. The same is true for antennas 207, 209, with one being operated at 400 kiloHertz and the other being 
operated at 2 megaHertz. However, the processors 215, 217. 219 can be programmed to provide any particular fre- 
quencies for the transmitting antennas. This will be used to good advantage as will be described below in connection 
with a calibration routine. 

In operation, a cornrnand signal is supplied to the FSELECT pin 351 to determine which frequency will be utilized 
55 for energizing a particular transmitting antenna. The FREQO register 327 and FREQ1 register 329 may be preloaded 
> ~ with two particular frequencies (such as 400 kiloHertz and 2 megaHertz). The binary signal applied to the FSELECT — - - — 
pin 351 determines the operation of multiplexer 353, which supplies the contents of either FREQO register 327 or 
FREQ1 register 329 of the input of phase accumulator 301 . Phase accumulator 301 accumulates a phase step on each 
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clock cyda The value of the phase step determines how many dock cycles are required far the phase accumulator to 
count two n radians, that Is, one cycle of th utput frequency. The output frequency is determined by the phase step 
multiplied by the frequency of the signal applied to the clock input pin 307 divided by 2 s2 . In practice, the phase accu- 
mulator 301 is cleared, then loaded with the output of multiplexer 353. Then, a predefined time interval is allowed to 
pass, during which the signal applied to clock input pin 307 steps the output of phase accumulator 301 through an incre- 
mentally increasing phase for the particular frequency. In other words, phase accumulator steps from 0° phase to 180° 
for a particular frequency. At any time, the output of phase accumulator 301 may be altered by a phase offset which is 
suppfied by phase register 331. Phase register 331 may be loaded in response to commands from processors 215. 
2 1 7. 21 9. The phase value is supplied as input to look-up table 303. which converts the output of the phase accumulator 
301 (and any desired offset) Into a digital bit stream which is representative of an analog signal. This digital bit stream 
is suppfied as an input to the 10-bit digital-to-analog converter 305 which also receives amplitude information from full 
scale control 359. The digital-to-analog converter 305 supplies an analog output with a particular frequency attribute, 
phase attribute, and amplitude attribute. For example, an output of 2 megaHertz, with 15° of phase, and a particular 
peak amplitude current may be provided as an input to a particular transmitting antenna. 

Figure 6 is a block diagram view of the digital signal processor 221 of Figure 3. In the preferred embodiment of the 
present invention, digital signal processor 221 comprises a DSP rrticrocomputer manufactured by Analog Devices of 
Norwood, Massachusetts, which is identified as Model No. ADSP-2101. This is a single^hp rnicrocomputer which is 
utilized for high-speed numeric processing applications. Its base architecture 379 is a fully compatible superset of the 
ADSP-2100 instruction set The base architecture includes three independent computational units: shifter 371, multi- 
plier/accurnulator 3? 3. and arithmetic and logic unit (ALU) 375. Program sequencer 369 supports a variety of opera- 
tions inducing conditional jumps, subroutine calls, and returns in a single cycle. Data address generator 367 includes 
two address generators. Digital signal processor 221 includes serial port 381 which includes two input channels: input 
channel 383, and input channel 385. Timer 387 provides timing signals for the data processing operation, and receives 
as an input a dock signal from divide-by circuit 239 (of Figure 3). External address bus 289 and external data bus 391 
allow digital communication between digital signal processor 221 and central processor 315 of Figure 6. Memory 393 
indudes program memory 395 and data memory 397. As is typical with digital signal processors, data memory 397 
defines at least two circular buffers associated with serial ports 383, 385. which are designed to receive asynchronous 
digital data, and store it indefinitely or for a predetermined time interval. The digital signal processor 221 receives digital 
inputs at channel inputs 383. 385 from an analogs-digital converter, such as is depicted in the drcuit of Figure 6. The 
receiving circuit of Figure 6 receives a current which is representative of the response of a particular receiving antenna 
21 1 , 213 to electromagnetic radiation propagating through the borehole. This electrical signal is processed through the 
circuit components of Figure 6. and is provided as an input to cBgital signal processor 221 . In the preferred embodiment 
of the present invention, receiving antenna 21 1 is identified with a particular input channel of digital processor 221 , 
while receiving antenna 213 is identified with the other input channel of digital signal processor 221. Central processor 
215 (of Figure 3) utilizes external address bus 389 and external data bus 391 to address a particular input channel and 
read digital data into central processor 21 5 tor processing. In the preferred embodiment of the present invention, digital 
signal processor 221 can sample data from receiving antennas 211. 213 at a very high sampling rate, which can be 
read periodically by central processor 21 5 which processes the data to determine the amplitude attenuation and phase 
shift of the electromagnetic signal which is propagated through the borehole. One particular routine for calculating 
amplitude attenuation and phase shift is set forth in greater detail herebeiow, in connedion with a discussion of the error 
cancellation feature of the present invention. In broad overview, central processor 215 can pull a selected amount of 
data from each channel of digital signal processor 221. and from that data calculate the amplitude attenuation and 
phase shift of the electromagnetic wave as it propagates through the wellbore and past receiving antenna 21 1 and 
receiving antenna 21 3. In the preferred embodiment of the present invention, an upper transmitter transmits an interro- 
gating electromagnetic signal of a particular frequency which propagates downward past receiving antennas 211,213. 
Then, a particular one of lower transmitting antennas 207, 209 propagate an interrogating electromagnetic signal 
upward. Measurements from receiving drcuit 231 are stored in the input channels of digital signal processor 221, and 
read by central processor 21 5 in a manner which allows for the calculation of amplitude attenuation and phase shift 

Another important feature of the present invention arises from the fact that a predse energizing current can be uti- 
lized to energize a particular one of transmitting antennas 203. 205. 207, 209. This will establish the frequency attribute, 
phase attribute, and amplitude attribute of the electromagnetic interrogating signal. Therefore, a single receiving 
antenna can be utilized to make the measurement of the electromagnetic interrogating signal as it passes through the 
wellbore. The amplitude and phase of that interrogating signal can be recorded in memory, and compared with values 
in memory for the energizing current This allows a single receiving antenna to be used to provide an accurate measure 
of amplitude attenuation between that particular receiving antenna and the particular transmitting antenna, and the 
phase shift of the interrogating signal between the transmitting antenna and the receiving antenna. -Of course.- the- 
amplitude attenuation and phase shift of the electromagnetic interrogating signal as it passes through the formation is 
indicative of the resistivity of the welfoore and surrounding formation. 
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Figures 7A. 7B, and 7C provide high level flowchart representations of logging operations performed in accordance 
with the preferred embodiment of the present invention. Figure 7A depicts logic steps which ar performed by central 
processor 21 5. Figure 7B represents operations controlled by processors 21 7, 21 9. Figure 7C depicts op rations con- 
trolled by digital signal processor 221 and central processor 21 5. Th transmission operations begin at block 401 . Proc- 

5 essor 215 performs a calibration operation upon receiving antennas 21 1, 213, as will be discussed in greater detail 
elsewhere in this application. After the calibration operations are performed central processor 215 instructs processor 
21 7 to energize transmitting antenna 203 with a 400 WloHertz current Then, in accordance with block 407, central proc- 
essor 21 5 instructs processor 2 1 9 to energize transmitting antenna 209 with a 400 WloHertz current Next central proc- 
essor 215 instructs processor 21 7 to energize transmitting antenna 205 with a 2 megaHertz current in accordance with 

it, block 409. Then, in occurrence with block 41 1 , central processor 21 5 instructs processor 219 to energize transmitting 
antenna 207 with a 2 megaHertz current The process stops at block 413. In actual practice, transmission operations 
will be performed continuously over predefined intervals. 

Figure 7B depicts the control operations performed by processors 217, 219 to cause numerically controlled oscil- 
lators 223. 225 to energize particular transmitters. The process begins at block 415. It continues at block 41 7. wherein 

is the processor 217 or 21 9 clears the registers in numerically controlled oscillators 223 or 225 by providing the appropri- 
ate instruction. Then, in accordance with block 419, processor 217 or 219 loads a predetermined value to the FREQO 
register and the FREQ1 register. These values determine the frequency of the energizing current which is supplied to 
a particular transmitting antenna. Inert in accordance with block 421, processor 217 or 219 loads a predetermined 
phase value to the phase register of numerically controlled oscillator 223 or 225. Processor 21 7 or 21 9 then provides a 

20 binary command to the FSELECT input pin of numerically controlled oscillator 223 or 225 to select a particular fre- 
quency of operation. Then, in accordance with block 425, a particular time interval is allowed to pass. This lime interval 
determines how many cycles of energizing current are applied to a particular transmitting antenna. The process ends 
at software block 427. Typically, each time processor 21 7 or 219 is instructed by central processor 215 to energize a 
particular transmitting antenna, the steps of Figure 76 are performed. 

25 Figure 7C depicts in flowchart for the reception operations. The process begins at block 429. The process contin- 
ues at block 431 , wherein the current within receiving antennas 21 1 , 213 are sampled by receiving circuit 231 . Then, in 
accordance with block 433, these samples are loaded to the appropriate input channels 283, 285 of digital signal proc- 
essor 221. In accordance with block 435, central processor 215 fetches selected samples from the memory buffers 
associated with the digital signal processor input channels. In accordance with block 437, optionally, samples may be 

30 modified to offset for error components due to "miscalibration" of the antenna, which will be described in greater detail 
Isewhere in this application. Next, in accordance with software block 439, the digital samples may be digitally f iter with 
either a low-pass digital filter, high-pass digital filter, or a bandpass digital filter. Alternatively, the samples can be aver- 
aged over predefined intervals to provide stability to the samples and eliminate the influence of spurious or erroneous 
samples. Next, in accordance with block 441 , the amplitude attenuation and phase shift are calculated, as is described 

35 elsewhere in this application. Finally, the process ends at block 443. 

Claims 

1 . A method of well-logging in a borehole in order to generate log data for a particular formation, comprising: 

40 

generating a model which maps, for a formation type similar to said particular formation, dielectric permitivity 
to anisotropic formation conductivity; 

making at least one measurement of said particular formation at a first interrogation frequency; 
making at least one measurement of said particular formation at a second interrogation frequency; 
45 utilizing said measurements and said model to calculate values of at least one of: 

(a) the horizontal and vertical conductivities of said particular formation; 

(b) the horizontal and vertical dielectric permitivity of said particular formation; and 

(c) a dip angle of said borehole. 

50 

2. A method of well-logging according to claim 1 . wherein said model is comprised of machine-readable data. 

3. A method of wed -logging according to claim 1 , wherein said steps of making at least one measurement comprise: 

65 making at least one measurement of said particular formation at one particular interrogation frequency; and 
making at least two measurements of said particular formation at a second particulaMnterrogation frequency. 

4. A method of well-logging according to claim 1 , wherein said step of utilizing said measurements and said model to 
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calculate comprises: 

wherein said first interrogation frequency comprises a frequency that is relatively low in comparison to said 
second interrogation frequency; 

utilizing said measurements and said mod I values to calculate at least one of th following: 

(a) the horizontal and vertical conductivities of said particular formation; 

(b) the horizontal and vertical dielectric perm'rtivity of said particular formation; and 
(cj a dip angle of said borehole. 

5. A method of well-logging according to claim 1, wherein said step of utilizing said measurements and said model 
value to calculate comprises: 

utilizing said measurements and said model values to calculate at least one of the following utilizing a four- 
dimension inversion: 

(a) the horizontal and vertical conductivities of sab particular formation; 

(b) the horizontal and vertical dielectric penmitrviry of said particular formation; and 

(c) a dip angle of said borehole. 

6. An apparatus for well-logging in a borehole in order to generate log data for a particular formation, comprising: 

a programmable controller; 

machine readable data which defines a model which maps, for a formation type similar to said particular for- 
mation, dielectric permitrvfty to anisotropic formation conductivity accessible by said programmable controller; 
a transmitter and receiver assembly for: 

(a) making at least one measurement of said particular formation at a first interrogation frequency; 

(b) making at least one measurement of said particular formation at a second interrogation frequency; 

program instructions for utilizing said measurements and said model to calculate values of at least one of: 

(a) the horizontal and vertical conductivities of said particular formation; 

(b) me horizontal and vertical dielectric permrtivity of said particular formation; and 

(c) a dip angle of said borehole. 

7. An apparatus for well-logging, according to claim 6, wherein said model comprises machine readable data. 

8. An apparatus for well-logging according to claim 6. wherein said transmitter and receiver assembly make at (east 
one measurement at one particular interrogation frequency and at least two measurements at a second particular 
interrogation frequency. 

9. An apparatus for wen-fogging according to daim 6, wherein said program instructions include program instructions 
for determination of particular ones of horizontal and vertical conductivities and horizontal and vertical dielectric 
penmitivrties through an inversion process. 

10. A method of welMogging in a borehole, comprising: 

providing a dip angle for said borehole; 

making at least one measurement of said particular formation at a first particular interrogation frequency; 
making at least one measurement of said particular formation at a second particular interrogation frequency; 
utilizing said measurements to calculate values of at least one of the following: 

(a) horizontal conductivity; 
ft)) vertical conductivity; 

— — (^horizontal dielectric permitivity; and - — - — — 

(d) vertical dielectric permftivfty. 
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1 1 . An apparatus for well-logging in a borehole, comprising: 

a programmable controller; 
a source for providing a dp angle for said borehol ; 
s a transmitter and receiver assembly for: 

(a) making a plurality of measurements of said particular formation at a first particular interrogation fre- 
quency; 

(b) making a plurality of measurements of said particular formation at a second particular interrogation fre- 
10 quency; 

program instructions for utilizing said measurements to calculate values of at least one of: 

(a) horizontal conductivity; 
is (b) vertical conductivity; 

(c) horizontal dielectric perrnitivity; and 

(d) vertical dielectric permitivHy. 

12. A method of well-logging in a borehole in order to generate log data for a particular formation, comprising: 

generating a model which maps, for a formation type similar to said particular formation, dielectric perrnitivity 
to anisotropic formation conductivity; 

making a plurality of measurements of said particular formation at a particular interrogation frequency; 
utilizing said measurements and said model to calculate values of: 

(a) horizontal conductivity; and 

(b) horizontal dielectric perrnitivity. 

13. A method of well-logging in a borehole, according to claim 12, wherein said step of utilizing said measurements to 
30 calculate includes performing a plurality of inversions to determine particular ones of horizontal conductivity; verti- 
cal conductivity, horizontal dielectric permrtrvity. and vertical dielectric permrtivity. 

14. A method of well-togging in a borehole, according to claim 12, wherein said plurality of inversions comprise a four- 
dimensional inversion performed as a series of two-dimensional inversions. 
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1 5. An apparatus for well-logging in a borehole in order to generate log data for a particular formation, comprising: 



a programmable controller; 

machine readable data which defines a model which maps, for a formation type similar to said particular for- 
40 mation, dielectric perrnitivity to anisotropic formation conductivity; 

a transmitter and receiver assembly for making a plurality of measurements of said particular formation at a 
particular interrogation frequency; 

program instructions for utilizing said measurements and said model to calculate values of: 

45 (a) horizontal conductivity; and 

(b) horizontal dielectric constant. 

16. An apparatus for welHogging in a borehole according to claim 15, wherein said program instructions for utilizing 
said measurements include program instructions for performing a plurality of inversions to determine at least one 

so of horizontal conductivity, vertical conductivity, horizontal dielectric perrnitivity. and vertical dielectric perrnitivity. 

17. An apparatus for welHogging in a borehole, according to claim 15, wherein said plurality of inversions comprise a 
four-dimensional inversion performed as a series of two-dimensional inversions. 
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